Possibility of Skyrmion Superconductivity in Doped Antiferromagnet K2Fe4Se5 
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Intercalated Fe chalcogenide K2Fe4Se5 family exhibits high Tc 30 K) superconductivity and 
spin-8 high Tn 560 K) antiferromagnet ism (AFM). We present a model Hamiltonian and suggest 
Skyrmion superconductivity. A doped electron creates an orbitally non degenerate S=7| state in 
a F4Se9 cluster and moves in a single correlated band, exchange coupled to a robust S=8 AFM 
order. Skyrmion, a topological excitation of 2 dimensional AFM order acquires an induced charge 
-2e through a quantum anomaly and becomes a Cooper pair. Superconductivity emerges for a range 
of doping. Fluctuating superconductivity, arising from preformed Cooper pairs (stable Skyrmions), 
around room temperatures is likely. There is hope for higher Tc's in other high spin Mott insulators. 
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Discovery of superconductivity in Fe based pnictides 
[U 12] and chalogenides [3] have created another wave of 
excitement in condensed matter physics. Even though Tc 
has not matched the scale of 100 K seen in cuprates, Fe 
systems continue to offer surprises. A recent surprise is 
the discovery [4 of superconductivity in intercalated Fe 
chalogenide family, where a robust antiferromagnetic or- 
der with a sublattice moment close to 13 /iB and a T^r ^ 
560 K seems to coexist with a high Tc superconductiv- 
ity of Tc ~ 30 K. From materials science point of view 
the vacancy order present in K2Fe4Se5 has enabled a rich 
physics. Ordered vacancies is a new route to novel states 
of quantum matter [5j. 

In the present article we provide a minimal model and 
identify a mechanism of electron pairing. We use some 
general considerations and present certain consequences 
that are compelling. We show that a single band model of 
strongly correlated electrons, in a square lattice, close to 
half filling exchange coupled to a square lattice of large 
spin-7^ core moments of Fe4 plaquettes captures the 
essential low energy physics. In our theory topological 
excitations of an ordered 2 dimensional Heisenberg anti- 
ferromagnet, namely Skyrmion [6], plays a fundamental 
role. Relevance of topological excitations [7], (FJ to lightly 
doped quantum antiferromangets began in 1987 with the 
advent of RVB theory of high Tc superconductivity. This 
study continues [TOl [11] . 

Briefly, kinetic energy frustration of doped electrons 
in the presence of a rigid AFM order is relieved by a 
Skyrmion which help delocalize 2 electrons, using its spin 
twist. Skyrmion, through a quantum anomaly, acquires 
an induces charge -2e and becomes a Cooper pair. For 
very low doping we get an ordered Skyrmion crystal. Be- 
yond a range of doping Skyrmion crystal quantum melts, 
leading to Skyrmion superconductivity. For even higher 
dopings ferromagnetism is likely to intervene, with a pos- 
sible island of p-wave superconductivity (figure 4). 

Following the discovery of superconductivity in FeSe 
[3], a new family, A^^Fe^ySe, with intercalated alkali atoms 
(A = K,Cs,Rb) was synthesized Stoichiometric 




FIG. 1: V5 a) vacancy (o) order and S=8 (=2 x 4) AFM 
order of Fe4 plaquettes in FeSe layer of K2Fe4Se5 . Only Fe 
atoms (•) are shown. Fe^^ is in S = 2, the maximal spin 
state, b) Se atoms above (red) and below (pink) Fe layers are 
shown. Dashed square in the middle marks a Fe4Se9 unit. 



K2Fe4Se5 has ordered Fe vacancies and y/b x peri- 
odicity. Experiments indicate [121 US] that K2Fe4Se5 is 
a Mott insulator. Neutron scattering [12] shows that the 
Fe4 plaquettes attain a large cluster moment and reaches 
a S=8 spin state. Further the plaquettes develop a (tt, tt) 
long range AFM order (figure 1) with a Tn ~ 560 K. 

Superconducting Tc is around 30 K, in various mem- 
bers of the family. Many experiments show coexistence of 
magnetism and superconductivity [Ql[13l[T5 . However, 
NMR and STM results [14 seem to challenge coexistence. 

There are some theoretical studies on this system |^ . 
We present a new approach and develop a minimal mi- 
croscopic model to understand magnetic and supercon- 
ducting properties. We first propose that Fe4 plaquette 
is a basic electronic building block. Electronic integrity of 
Fe4 plaquette is revealed by its robust total spin-8 state 
arising from strongly exchange coupled four Fe atoms 
that carry spin-2 each. The S=8 plaquette spins have 
a (7r,7r) AFM order in the 2D square lattice of plaque- 
ttes. An important question is where do doped carriers 
go ? This question is similar to Zhang-Rice singlet for- 
mation in cuprates. In view of available experiments we 
will focus on electron doping. 

Simple quantum chemical and phenomenological con- 
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FIG. 2: a) Sciiematic energy levels of 16 singly occupied and 
spin polarized plaquette orbitals of the Fe4Se9 cluster, b) 
The strongly hybridized core spin- 7 1 moment and the valence 
spin- 1 moment distributed in the Fe4 plaquette is schemati- 
cally shown. 



siderations suggest that a doped electron enters a 
Fe4 plaquette and create an orbitally non degenerate 
spin-7^ , hexadectic state in the Fe4 plaquette [16] . 
Briefly, Fe atoms of Fe4 plaquette are tetrahedrally co- 
ordinated by Se atoms and form a Fe4Se9 unit. In view 
of strong hybridization of 3d orbitals of Fe and 4s/4p or- 
bitals of Se atoms, the 16 electrons that makeup S=8 spin 
state of Fe4 plaquette reside in 16 plaquette orbitals in a 
strongly exchange coupled parallel spin state (figure lb). 
A doped electron enters the top most singly occupied 
non degenerate orbital. We call this as valence orbital of 
Fe4 plaquette that carries a spin-^ valence moment. 

Anticipating addition of an electron (doping) we conve- 
niently separate S=8 spin of Fe4 into S=7^ core spin mo- 
ment Si, ferromagnetically coupled by a plaquette Hund 
coupling Jh to a S=^ valence moment (figure lb). 
Neighboring Fe4 plaquette spins are antiferromagneti- 
cally coupled by an effective superexchange J. Thus the 
spin Hamiltonian of the Mott insulator is: 

i^AFM = -JH^S,-S, + J^(S,+S,)-(S, + S) (1) 



{ij) 



Here = cJ^aQ,/^^^, is the valence spin-^ operator writ- 
ten interms of electron operators (cj^, c^^) of the Fe4 pla- 
quette valence orbital. For simplicity we ignore the weak 
c-axis exchange coupling and any anisotropic spin inter- 
action terms. 

In the limit Jh ^ J, equation 1 becomes a Spin-8 
Heisenberg AFM Hamiltonian on a square lattice. The 
effective exchange coupling J can be estimated from the 
acoustic spin wave band width from inelastic neutron 
scattering results [12 . The estimate is JS = 8J 17 
meV. Location of the first optic spin wave branch give a 
value for Jh ^ 110 meV. 

The Hamiltonian of our electron doped system is a sin- 
gle band model of plaquette valence spin electrons close 
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FIG. 3: Spin configuration of a Skyrmion. A circle in the 
middle has only xy component of sublattice spins and defines 
an XY- vortex. An anti Skyrmion will have an antivortex. 
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with a local constraint (^a^ia ^ (^^^ electron dop- 
ing). Here t is the inter plaquette hopping matrix ele- 
ment. From known band structure results, after folding 
the BZ due to vacancy ordering, we estimate t ~ 50 meV. 
The band width W ^ 400 meV is large compared to the 
charging energy Ueff ^ ^ ~ 500 meV of the (Fe2Se5)^+ 
plaquette. Here R is the size of the Fe4 plaquette. Our 
doping fraction x is with reference to the Fe4 plaque- 
tte orbital band. Thus in the formula K2+yFe46'e5, our 

doping X = ^. 

The low energy effective Hamiltonian of the interacting 
core spin Heisenberg antiferromagnet is a 2 + 1 dimen- 
sional 0(3) nonlinear sigma model with action: 



dx dy dT[{dij,n{r)Y 



(3) 



Here /i = x,y; Vg is the spin wave velocity and n(r) is the 
normalized, n(r) • n(r) = 1, sublattice magnetization. 

The above system supports small amplitude gapless 
spin wave (Goldstone mode) excitations and finite energy 
topological excitations namely Skyrmions. Skyrmions 
are analogue of a pair of domain walls of ID Heisenberg 
antiferromagnet. In ID system a single domain wall is a 
TT— twist that interpolates up and down spin vacua in a 
continuous fashion. In the process it effectively decouples 
a single spin (a spinon localized around the domain wall) 
from the antiferromagnetic mean field. In the language 
of spin density wave (similar to polyacetylene) it creates 
a midgap state containing an unpaired spin. 

A Skyrmion interpolates the down spin vacuum in its 
central region and up spin vacuum at the boundary, in 
a radially symmetric fashion (figure 3). In the process it 
effectively decouples two spin- 1 spinons [8 . This can be 
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seen in several ways [7'-^. A simple way is in the SDW 
band picture, as studied by John and collaborators [9 . A 
skyrmion or 2 half skyrmions (2 merons) configuration 
traps two unpaired spins in two mid gap states. 

Skyrmion has topological properties. It carries chiral- 
ity or Pontryagin index Q (= 0, ± 1, ± 2, ...). Its sign 
measures the sense, and magnitude number of windings, 
when we map order parameter n(r) field from the 2D 
plane onto S^, the surface of an ordinary sphere: 

Q = J dx dy n- {dxU x dyii) (4) 

In the language of homotopy theory, 772 (S^) = Z. 
Skyrmion has a zero mode associated with a residual 
U(l) symmetry about the z-axis of the spin configura- 
tion. Quantized zero mode represents spin current fiuc- 
tuation [17 and quantum Skyrmion is a chiral doublet. 
This is a 2D analogue of Affleck's dyon[T8] in ID Heisen- 
berg antiferromagnets. 

Valence electrons living in a single band in our model 
are exchange coupled to an antiferromagnetic mean field 
provided by the core spin moments from Fe4 plaqauettes. 
An electron in a rigid AFM background can not delocal- 
ize. We have kinetic energy frustration. Skyrmion (Q = 
± 1) relieves the frustration by delocalizing 2 electrons 
using spin current (spin twist) under its wings. 

Taking care of the local constraints and exchange in- 
teractions in a mean field fashion we get the mean field 
Hamiltonian of the valence electrons system as: 

i^mF ^ ^ij4a^jcr + ^.C - J^S ^ 11^ • cj^aa^^/? (5) 

ij i 

Here iij is a renormalized interplaquette hopping ma- 
trix element, (we choose tt fiux mean field solution); 
Jh = Jh + 4 J. The (tt, tt) antiferromagnetic order in- 
troduces a gap in the two Dirac nodes. To calculate in- 
duced electric charge, in the presence of a slowly varying 
order parameter n(r) and U(l) gauge field A^, we use 
the language of 2D Dirac fermions and follow Wigmann 
and Abanov[7]. The effective action of the valence elec- 
tron quasi particles coupled to slowly varying background 
AFM field is: 

S ^ J d^xip [-i-ff,{d^ -^Al^)^ima- n(r)] (6) 

Here /i = x,y and ip = iip^j^, and 7^,7^,7^ = 
a^^a^ , —cr^. We have rescaled the Fermi velocity of elec- 
tron to 1 for brevity. Further m ^ 6'Jh is a mass gap 
provided by the strong exchange field. Integrating out 
the fermions and in the large mass approximation one 
gets ^: 

S = S^ + 2iA1Jl (7) 

Here Jj(r) = ^e^^'^n - djyW x d\w is the topological 
spin current density, whose time component Jq{v) = 



• {dxn X dyu) is the Skyrmion density (equation 4). 
The real part of the induced interaction renormalizes the 
NLS action (equation 3). The imaginary part is the key 
U(l) gauge field - topological spin current Jj cou- 
pling. The time component of the above, 2zAq Jq^ is the 
coupling of dopant chemical potential to Skyrmion den- 
sity (equation 4). Hence we get an induced U(l) (electric) 
charge Qskyrmion given by the expression [7]: 

Qskyrmion = ± 2e \Q\ (8) 

where Q is the Pontryagin index or Skyrmion number. 

To calculate the spin of the - 2e Skyrmion, we start 
with an even number of electrons in an ordered antifer- 
romagnet and create one Skyrmion. This generates two 
singly occupied mid gap states. Addition of two electrons 
to these midgap states completely fill the two states and 
create spin zero midgap states. Thus nominally a 2D 
Skyrmion carrying charge - 2e is a spin singlet. 

2e Skyrmion, a charged chiral doublet, as discussed 
by Lee [11 for the case of a spin liquid, represents a 
= ± 2 angular momentum, d ± id state of a Cooper 
pair. Skyrmions that are real combination of the above 
doublet will lead to a real d-wave superconducting state. 
Otherwise we will get a d ± id superconducting state that 
breaks PT symmetry. This will be a chiral superconduc- 
tor d la Laughhn. 

In an ideal 0(3) non linear sigma model in 2D, 
Skyrmion energy is finite and size independent. In 
practice, coulomb repulsion among 2 electrons, spin 
anisotropy and interchain coupling determine the size. 
For our case, the mean size of Skyrmion is determined 
by the inter Skyrmion distance. Cooper pair size or co- 
herence length is selfconsistantly limited by mean inter 
Cooper pair distance. 

Skyrmion delocalization and kinetic energy gain is de- 
termined by their mutual overlap, proportional to x. As 
a first approximation, Skyrmion band width is a fraction 
X of the fermi energy Wx of free dopants. That is, the 
Skyrmion band width Wskyrmion ~ Wx^. Since W ^ 400 
meV, Wx^ ~ 4 meV, for a doping of x ~ 0.1. As -2e 
Skyrmions are Bosons this is also Bose condensation en- 
ergy scale. This energy is in the right ball park of the 
observed Tc ~ 30 K. Our superconductor is best viewed 
as a 2D Kosterlitz-Thouless condensed state of interact- 
ing bosons. Weak Josephson coupling between layers will 
make it a 3D superconductor. 

Skyrmions are stable, in principle upto Neel temper- 
ature, as long as AFM order survives. In this sense 
2e Skyrmions are preformed Cooper pairs in the back- 
ground of an AFM order. As a consequence we could see 
Cooper pairs well above Tc, through for example anoma- 
lous Nernst effect. We predict that K2Fe4Se5 family of- 
fers a wonderful and some what unique opportunity to 
explore and utilise ffuctuating superconductivity around 
room temperature scales. 
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FIG. 4: Schematic phase diagram for doped K2Fe4Se5 family. 



In our theory the large spin AFM order and super- 
conductivity coexist for a range of doping. Dopant de- 
localization will eventually destroy the long range AFM 
order. In the low doping regime unscreened coulomb in- 
teraction will localize Skyrmions resulting in Skyrmion 
crystal or Skyrmion glass. Berciu and John [9^ have dis- 
cussed this for cuprates. A universal resistivity maxima 
seen in K2Fe4Se5 superconductor could indicate onset of 
Skyrmion crystallization frustrated by growing quantum 
fluctuations (dashed line in figure 4 inside AFM dome). 

It is known that large spin doped antiferromagnets 
leads to ferromagnetism through double exchange^ for 
example in doped manganites, a spin-2 Mott antiferro- 
magnet. As shown by de Gennes, small doping leads to 
spin twist (canting) , before ferromagnetic order sets in at 
higher dopings. In a 2 dimensional system such as ours, 
the system naturally organize spin twists in the form of 
topologically stable Skyrmion eigen excitations and help 
delocalization. What is exciting is that superconductivity 
is an alternative to ferromagnetism in 2 dimensions. So 
we have a new avenue for high Tc search, namely high 
spin antiferromagnets in lower dimensions. 

It also follows from double exchange physics that heav- 
ily doped K2Fe4Se5 will eventually become a ferromag- 
netic metal. If we take clue from Sr2Ru04, another 2D 
even electron (Ru^+ = Ad^) system, Hund coupling and 
a growing ferromagnetic fluctuation could stabilise an is- 
land of spin triplet p-wave superconductivity. 

A robust AFM magnetic order of a large spin Mott 
antiferromagnet offers exciting possibilities, as shown 
schematically in figure 4. It will be useful to study doped 
K2Fe4Se5 family experimentally with controlled doping. 
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